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Abstract. We investigate a periodic quantum graph in form of a square lattice with 
£Nj . a general self-adjoint coupling at the vertices. We analyze the spectrum, in particular, 

its high-energy behaviour. Depending on the coupling type, bands and gaps have 
different asymptotics. Bands may be flat even if the edges are coupled, and non-flat 
T^j- . band widths may behave like 0(n J ), j = 1, 0, — 1, —2, —3, as the band index n — >• oo. 

t-H ' The gaps may be of asymptotically constant width or linearly growing with the latter 

\Q ■ case being generic. 

o : 
o : 

' 1. Introduction 

Origins of the quantum graph concept can be traced back to Linus Pauling's 
considerations about the structure of organic molecules in 1930's, however, it was 
roughly the last two decades when the subject became very popular and its true richness 
has been revealed. The reason for that was not only the fact that quantum graphs 
represent a suitable model of numerous systems prepared by solid state physicists. 
Equally important was the inherent structure of such models which allowed us to 
study effects uncommon in the "usual" quantum mechanics coming from the nontrivial 
topological structure of graph structures as well as from the fact that they mix features 
corresponding to different dimensionalities. The bibliography concerning quantum 
graphs is nowadays extensive indeed and we restrict ourselves to quoting the proceedings 
volume of a recent topical programme at Isaac Newton Institute and references therein 
|AGA08| . 

One often studied class are periodic graphs. Their spectrum has predictably a band 
structure, however, in distinction to the usual Schrodinger operators they can exhibit 
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under particular geometric conditions also infinitely degenerate eigenvalues manifesting 
invalidity of the unique continuation principle |Ku05] . The aim of the present paper 
is to investigate the spectrum in a particular example of periodic graphs, namely 
two-dimensional square lattices with a general self-adjoint coupling at the vertices. 
This generalizes earlier work in which lattices with 5 and 5' s couplings were studied 
[Ex96llEG9^ . 

A motivation for the present extension comes from different sources. First of all, 
the general vertex coupling became more interesting after several recent results - 
cf. |KZ014 IEP09t ICET10] and references therein — showing how it can be approximated 
by graphs or network systems with suitably scaled potentials illustrating thus that it is 
not just a mathematical object but also something which can be, in an approximative 
sense at least, realized physically. Secondly, already the mentioned examples of 5 and 5' s 
couplings demonstrated that different couplings yield different asymptotical behaviour 
of spectral bands and gaps at high energies. 

These asymptotics can have interesting dynamical consequences. Recall that in 
the one-dimensional analogue of the present problem, in the generalized Kronig-Penney 
model, there are three types of asymptotic behaviour [EG99, CS04] and that the one 
corresponding to the 5' interaction for which spectral gaps are dominating exhibit 
absence of transport in the Wannier-Stark situation when an electric field is applied 
[AEL941 lEx95l IMS951 KDE98] : we stress that the solution of the Wannier-Stark problem 
in the other two cases mentioned is still an open question. With these facts in mind 
it is natural to ask how many types of asymptotic behaviour a two dimensional square 
lattice can show and what they look like. Our approach to this problem is based on the 
use of the so-called ST-form of boundary condition introduced in |CET10j . This will 



allow us to express the band and gap widths by rather simple expressions involving the 
coupling parameters, from which it is easy to determine how the spectral behaviour is 
governed by the particular boundary conditions. 

We will show that the high-energy behaviour has more types than in the one- 
dimensional situation. In fact, the problem has sixteen parameters and offers a zoology 
of solutions. The scope of this issue does not allow us to present a complete classification 
but we will list all the "generic" cases with respect to the rank of the matrix B in the 
condition (CQ) below and single out cases of particular interest. The two-dimensional 
character of the problem has two main consequences. The first one is a possible 
occurrence of flat bands, or infinitely degenerate eigenvalues; this is connected with 
the invalidity of the principle of unique continuation on graphs of nontrivial topology 
mentioned above. The second remarkable feature is the existence of couplings for which 
the spectral bands are shrinking as n^ 1 , n~ 2 , or n~ 3 with respect to the band index 
n. This effect again has no analogue in the one- dimensional situation. The most 
"generic" situation, however, is the one known from [Ex96, EG96j when the bands are 
asymptotically of constant widths and gaps are linearly growing. 
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2. Preliminaries and main result 

2.1. Square lattice with a general vertex coupling 

As we have said we will consider a square lattice graph. To be concrete, the vertices of 
r are {(ma,na) : m, n G Z} for a fixed a > and the edges are segments of length a 
connecting points differing by one in one of the two indices. The state Hilbert space is 
the orthogonal sum of the L 2 spaces on the edges and the Hamiltonian acts as — -P4 on 
each of them, with the domain consisting of the corresponding W 2,2 functions. 

It is well known that in order to get a self-adjoint operator one has to impose 
boundary conditions at graph vertices which couple the vectors \I/(0) and ^/'(O) of the 
boundary values — we choose the variables at all the adjacent edges so that they start 
at the vertex. The standard form of these conditions is 



A$(0) + BV'(0) = 0, (1) 

where A, B are matrices such that (A\B) has maximum rank and AB* is self-adjoint 
[KS99J. Being interested in the periodic situation we naturally assume that the matrices 
A, B are the same at each vertex of the lattice. A drawback of the conditions (CQ) is 
that the matrix pair is not unique. There are various ways to mend this problem. One 
is to rewrite (CD) in the form 

(U-I)V(0)+i(U + I)V'(0) = 0, (2) 

where U is a unitary matrix. In the quantum graph context these conditions have been 
proposed in |HaOQ[ [KS00], however, they were known earlier in the general theory of 
boundary forms [GG91J. It is clear that a distinguished role is played by subspaces of 
the boundary space values referring to eigenvalues =Fl of U. An alternative way |Ku04] 
to write the conditions is by means of the corresponding orthogonal projection P and 
its complement Q := I — P: there is a self-adjoint operator L in QC n such that 

Ptf(0) = 0, Qtt'(0) + LQ^(0) = 0. (3) 

Here we are going to use yet another version introduced in |CET10] as the ST -form, 



T o ) no) = ( -t> A ) * (0) (4) 

using matrices T G C m ' n_m and a self-adjoint S G C m ' m , where n is the dimension of the 
boundary value space (which will be four in our case) and m = n — dimP G {0, . . . , n}. 



2.2. The main result 

Let Hs,t is the quantum graph Hamiltonian described above. Our results about its 
spectrum can be summarized as follows: 



High- energy asymptotics of the spectrum of a square lattice 



4 



Theorem 2.1. (a) The spectrum of Hs.t consists of absolutely continuous spectral bands 
and infinitely degenerate eigenvalues. Its negative part consists of at most four bands. 

(b) The high-energy asymptotic behaviour of spectral bands and gaps as a function of 
the band index n includes the following classes: 

• flat bands, i. e. infinitely degenerate point spectrum, 

• bands behaving like 0{inP), j = 1, 0, —1, —2, — 3 ; as n — » oo, 

• gaps behaving like 0{vP), j = 1,0, as n — > oo. 

Depending on the vertex coupling (CQ) the high- energy asymptotics of the spectrum may 
be a combination of the above listed types. 

The rest of the paper is devoted to demonstrating of these claims. We will do that 
by analyzing spectrum of the fiber operator coming from the Bloch-Floquet analysis, 
discussing subsequently situations corresponding to different values of rank m of the 
matrix B in the boundary conditions (TjQ). While this procedure serves best our aim, 
it is not very illustrative from of the point of view of particular type of edge coupling. 
Apart from the trivial case m = when the lattice decomposes into separate edges 
with Dirichlet conditions, each of the other values of m cover several subcases with very 
different couplings and spectral behaviours. They would deserve a separate discussion 
which we cannot present here due to volume restrictions and we postpone it to another 
publication; we limit ourselves to several general statements: 

• The generic situation corresponds to m = 4 with all the edges coupled and spectral 
gaps growing linearly with the band index, 

• each case m = 1, 2, 3, 4 covers a situations with flat bands corresponding to lattice 
decoupling into separated edges, or pairs of edges, 

• the lattice can separate into "one- dimensional" subsets describing generalized 
Kronig-Penney models on lines or zigzag curves, or to "combs", 

• from the spectral point of view the case m = 3 is the richest, including situations 
with a powerlike shrinking of spectral bands that occurs for the graph decomposed 
into "combs". 

3. Bloch-Floquet analysis 

Since our graph is a-periodic w.r.t. shifts in both directions, we are able to employ the 
Bloch-Floquet decomposition. The elementary cell is depicted in Fig. [TJ together with 
the notation of the wave function components on the edges. 

The analysis follows the same pattern as in [Ex96, EG96] but the general vertex 
coupling makes it substantially more complicated. The fiber operator corresponding to 
fixed values of Floquet parameters (quasimomentum components) has a purely discrete 
spectrum and the number of its negative eigenvalues is at most four. The last claim 
follows from general principles [We80, Sec. 8.3] and from a comparison to the square 
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Figure 1. 4 periodic two-dimensional network 



graph Hamiltonian with Kirchhoff coupling at the vertices the spectrum of which equals 
by [Ex96j it is enough to notice that the Floquet component of this operator and 
that of a general H$,t have a common symmetric restriction with deficiency indices 
(n, n), n < 4. Note also that the bands may overlap in general; an example can be 
constructed using boundary conditions which separate motion in the two directions 
leading to two families of generalized Kronig-Penney models. 

Our main interest concerns the positive part of the spectrum, and as usual we are 
thus going to investigate solutions of the corresponding stationary Schrodinger equation 
with energy E = k 2 , k > 0. It is obvious that they are at each edge linear combinations 
of the functions e lkx and e~ lkx , specifically 



■02 (z) 



C+ e ikx + Cfe 

'2 

Dfe ikx + Dfe 
D+e ikx + D^e~ ikx , 



-ikx 



x G [-a/2,0] 
( e ikx + C^e~ ikx , x G [0, a/2] 

x G [-a/2,0] 
x G [0, a/2] 



-ikx 



(5) 



By assumption, they have to satisfy the boundary conditions at the vertex, i.e. it holds 



.4 



/ ^i(O) \ 

02(0) 

¥>i(0) 

V <pM J 



+ B 



( -V4(o) \ 
V4(o) 

-^i(0) 

V ^ 2 (o) J 



0. 



(6) 



In addition to that, of course, they have to satisfy Bloch-Floquet conditions, 

Ma/2) = e iei ^(-a/2) p 2 (a/2) = e"V(-a/2) , 

^(a/2) = e i6 ^[(-a/2) ^(a/2) = e w ^[(-a/2) , 



(7) 



for fixed values of the quasimomentum components 61,62 G (—71,71]. Substituting (jHJ) 
into (J?D allows us to express the variables and Df in terms of Cf and Df, 



(8) 
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Using these relations we eliminate and Df from (jSJ), after that we substitute (151) 
into ([6]). Simple manipulations then yield the following condition, 



[(AM + \kBN)D] 



( ct \ 

\ d t J 



(9) 



where D := diag (e^ 6l ~ ak \ e ^ 6l+ak \ e W 2 - ak \ e ^ d2+ak) ) and the matrices M, JV are 
given by 



M :-- 



\ 



e i(0l-afc) 






/ _Q-^(Si-ak) 



N 



C2 



i(6»i-afe) 



V 



o 
o 



Q2 




e -5(^i+afc) 
_Q^(8i+ak) 






e 2 i 
e^ ! 






9 2 -afc) 
9 2 -afc) 



e 2 






i(0 2 +afc) 



02 



Ue 2 +ak) 






_ e -^(d2-ak) e -^(92+ak) 

Q^(02-ak) _Q^(d2+ak) 



J 



It follows from ()8]) that the functions ([5]) correspond to a nonzero solution iff 
(Ci,Ci,D^,Di^j is a nonzero vector. Consequently, a number A; 2 belongs to the 
spectrum of the Hamiltonian if and only if (Q has a non-trivial solution for some pair 
($i)#2), in other words, if there are 9i, 82 such that det [(AM + ikBN)D] = which 
can be simplified to 

det (AM + ikBN) = . (10) 

Our aim is to analyze the spectral asymptotics in dependence on the coupling type. 
Since four edges join at each vertex, the problem has 16 real parameters. We will take 
them into account through the ST form (jlj) of the boundary conditions, i.e. we set 





where m := rank(5). Obviously, each value of m has to be discussed separately. 



4. The case of m = 0, or Dirichlet decoupled edges 

Consider first the simplest situation when B = 0, A = —I. The ST-form of boundary 
conditions is then obviously invariant with respect to the edge labelling and the spectral 
condition f fTUj) acquires the form det(— M) = 0. Since det(— M) = det(M) = —4 sin 2 ak, 
this requires sinak = 0, hence the spectrum consists of infinitely degenerate eigenvalues, 



a(H) 



n7r\ • 
a ) 



n e N 
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5. The case of m = 1 

The admissible couplings form a seven-parameter family corresponding to the choice 

is \ 
-^100 
10 
\ -h 1 ) 

with s G M. and tj G C. Note that while in general the boundary conditions in the 
ST-form are not invariant with respect to the edge renumbering, we can choose without 
loss of generality the "privileged one" corresponding to the first component. A direct 
calculation of the determinant in (|T0l) leads to the spectral condition 

-4A;sinafc [(1 + |ti| 2 + |t 2 | 2 + |t 3 | 2 )cosa£; - 2$t(t 1 e Wl +ht 3 e W2 )] - As sin 2 ak = 0. 

We see that the condition is solved again by — for n G N. In addition to that, there 
are solutions coming from the equation 

(1 + |ti| 2 + |t 2 | 2 + \h?) cos ak - 23tf(W 01 + T 2 t 3 e ld2 ) = ~s sinak. (11) 

k 

Recall that k 2 > is in the spectrum if there are 61,62 G (— 7r,7r] which together with 
k satisfy (ITT]) . It is convenient to rewrite the range of the expression 3t(tie l61 + t 2 t 3 e 2 ) 
using a single parameter as follows: 

Observation 5.1. It holds 

{K(£ ie iei +r 2 t 3 e W2 ) I 6 l ,6 2 G (-7r,7r]} ={(|t x | + \t 2 \ ■ |t 3 |)cos0| $ G (-7r,7r]} . 

The spectral condition (ITTj) then yields the following requirement: 

(1 + |ti| 2 + \t 2 \ 2 + |t 3 | 2 )cosaA; -2(|ti| + \t 2 \ ■ |t 3 |) cos t9 = -~s sinafc. (12) 

This conditions has various types of solutions in dependence on the parameter values. 
Before discussing them, let us mention another useful fact. 

Observation 5.2. It holds 2{\t\\ + |t2| ■ 1^1) < 1 + |ti| 2 + | ^2 1 2 + \t 3 \ 2 and the equality 
occurs if and only if \t\\ = 1 A ^1 = l^l- 

5.1. Point spectrum 

If the expression \t\\ + \t 2 \ ■ \t 3 \ vanishes, i.e. if t\ = and t 2 = V t 3 = (we may 
suppose without loss of generality that t 3 = 0), then (ITT]) becomes 

cotgafc = -i r ^— . (13) 



/ 1 



B 



tx 














.4 
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The rhs being 0{k x ) as k — » oo the solutions are obviously close to the numbers 
(— \ +n)- at high energies. Writing them as 

k= {-\ + n )l + s - 

we get cotgak = a5 + 0(5 3 ) and t = — + 0{n~ 2 ). A substitution into f|T3|) and a few 
manipulations give 5 — ^ ■ 1+ ~^| 2 + O {n~ 2 ), hence this solution represents the spectral 
points k 2 = k 2 behaving like 



k 2 



1 \ 7T 

hn - 

2 ) a 



2 

a 1 + 1*2 



for n — > oo. These solutions are independent of the quasimomentum and give rise to 
flat spectral bands; it is not difficult to see that the corresponding eigenfunctions can 
be chosen compactly supported. 

Let us remark that the boundary conditions studied above mean that the 
Hamiltonian decouples into a countable direct sum of operators supported on two edges 
of the graph, or on individual edges if T = and s — 0, and therefore it is not surprising 
that the spectrum is pure point. 

5.2. Linearly growing spectral bands and gaps 

If | ii| + |t2 1 • | £3 1 7^ 0, we can divide by it and express thus cos^ from f fT2|) . This yields 
the spectral condition in the form 



1 + 


*i| 


2 + 


t-2 


to 




h 


2 


2( 




| + 


h 







, 1 s 

COS Ok + — • —r. : ; : : rr Sill dk 

k 2 ii + i 2 ■ i 3 



<1. (14) 



By Observation 15.21 the coefficient of cos ak at the Ihs cannot be smaller than one. 

In the rest of this section, we focus on the case |ii| 7^ 1 V \t-z\ 7^ \t$\\ the 
remaining situation will be treated in Section 15.31 Since the coefficient of sin ak is 
0(k^ 1 ), it is evident that dHj) can be asymptotically satisfied only away from the points 
where | cos ak\ = 1, in other words, the spectral bands are neigbourhoods of the points 
[(— I + n) ^ ] . Let us set k = (— ~ + n) - + d and find the range of d. We see that 
cos ak = sin ad and -r of the rt-th band solution is thus p4p can be rewritten as 



Hence ad < A := arcsin ^n^p^r^j^p i n °te that A G (0, |) in view of the assumption 
and Observation 15.21 Then the spectral bands behave like 



2(|*i 


| + 


to 




*s|) 


1 + 


*i| 


2 + 


*2 


2 


+ 


h 


2 
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which can be rewritten as 



I \ 2 7r 2 mr A m , . / 1 \ n~ ;/t , 



2 2 
7T 



in the high-energy limit, n — > oo. In other words, both bands and gaps are 
asymptotically linearly growing with the band index. 



5.3. Asymptotically constant spectral gaps 

If = 1 and \t2\ = \ts\ ytz 0, the spectral condition (Il4p has to be treated differently 
being now of the form 



cos ak + — 



k 2(l + |t 2 | 2 ' 



■ sin ak 



< 1 . 



(15) 



Let us suppose that s ^ — the case s = is special and will be discussed below 
in Section [531 Since | cosafc| < 1 and the coefficient at sin ak is small in modulus for 
large values of k, we see that the condition (Tl5l) is violated only in a small one-sided 
neighbourhood of the points where | cosafc| = 1. To describe the corresponding gaps, 
we set k = — + 5, then we have 



cos ak 



[aSf 



+ O (5 4 ) 



smak = (-l) n • a5 + O (S 3 ) 
1 a m ( 5 

+ o 



k nn 

Substituting from here into the negated condition ( fl5|) we obtain the gap condition, 



1 - 



(aS)* <r 



+ 



5 + 0(5 4 )+0[-)+0 /iy 



2 ' n7r2(l + \t 2 \ 2 j 
which is, in dependence of the sign of s, solved by 

Consequently, the gap boundaries are 



> 1 



5 e 



^Ly + Oin' 1 ) and 



/nn\ 2 I s 
\~a) + a ' 1 + \t 2 \ 2 



+ 0(n- 1 ) 



in other words, the gap widths are asymptotically constant and spectral bands grow 
linearly w.r.t. the band number n. Note that this case includes lattices with a nontrivial 
5 coupling discussed in |Ex96[ IEG96] . 
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5.4. No gaps, spectrum on the nonnegative halfline 

It remains to discuss the case with |ti| = 1, ^1 = l^l 7^ and s = 0, when the 
spectral condition ( TT5l) simplifies to |cosafc| < 1 which is obviously satisfied for all 
k > 0. Moreover, one checks directly that G cr(H), and putting k = in we find 
a(H) n (-00, 0) = 0, hence a(H) = [0, +00). Referring again to |Ex961 IEG96] we note 
that this includes the case of a lattice with Kirchhoff coupling. 



6. The case of m = 2 



The reader has noted already that the boundary conditions in the ST-form are not 
invariant with respect to the edge renumbering. Neglecting trivial lattice replacement 
corresponding to rotations and mirror images, we must distinguish two situations here: 

(i) Linearly independent columns of B are associated with parallel edges; without loss 
of generality we may suppose they are the "horizontal" ones. Then we apply the 
conditions in ST-form directly, 



B 



( 1 % 

1 t 21 



y 



tl2 

t 22 





\ 



/ 



A 



( 



\ 



sn S12 \ 

s 2 2 

-t^ 1 

-*12 -<22 1 ) 



where s u , S22 are real and the other matrix entries are complex. 
Linearly independent columns of B can correspond also to mutually orthogonal 
edges, say, the left "horizontal" and the lower "vertical". Then we use the conditions 
in a permuted form, with second and third row of \I/(0) and ^/'(O) interchanged. 
Since it is convenient to keep the entries order in these vectors, we interchange 
instead the second and the third column of the matrices A, B, 



B 



/ 1 til *12 \ 

£21 1 t 22 







A 



I s n \ 

S12 s 22 

-in 1 

~tl2 -taa 1 



However, since the spectral analysis of situations (i) and (ii) can be done in a very 
similar way, and moreover, also the structure of the spectral bands is essentially the 
same, we perform the analysis of the case (i) only. 

The determinant in (fTOj) leads to the spectral condition which can be written as 



V 2 ■ k 2 + V 1 ■ k + V = 
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where V 2 , V\ and V are expressions depending on ak, on the entries of S, T, and on the 
quasimomentum components 9i,9 2 . By a direct computation we get 

V 2 = -4cos 2 ak [\t u \ 2 + \t 22 \ 2 + \t 12 \ 2 + \t 21 \ 2 ] + 4sin 2 aA; [l + \tnt22 - t 12 t 21 \ 2 ] 
+ 8cosak [-^ ((t u t^ + t 12 t^)e Wl ) + 3? ((t 22 t^ + t^t 12 )e W2 )] 
+ 83? [tnt^e^-^] + 83? [*i 2 i^ 1+ ^] , 

Vi = 4sina/c [-cosaA; (s n (l + |t 2 i| 2 + \t 22 \ 2 ) + s 22 (l + \t u \ 2 + |t 12 | 2 ) 
-23? (512(^21+^22))) 

-23? (s 12 e i91 ) + 23? ((s u W22 + s 22 hit 12 - s 12 ~hit 22 - Si 2 t^2i)e ie2 )] , 
V = -4 sin 2 ak ■ det S . 

6.1. Linearly growing bands and gaps, or absence of gaps 

If we divide the above spectral condition by k 2 , we can write it in the asymptotic form 
V 2 {ak) = Oik' 1 ), explicitly 

- Acos 2 ak [\t u \ 2 + \t 22 \ 2 + \t 12 \ 2 + \t 21 \ 2 ] +Asin 2 ak [l + \t u t 22 - ti 2 t 2 i| 2 ] 
+ 8cosaA; [-3? ((t n t^ + t 12 t^)e Wl ) + 3? ((t 22 t^ + t^t 12 )e i92 )} 

+ 83? [h^ 9 ^] + 83? [t l2 hij {ei+92) ] = O (AT 1 ) , 

from which it is possible to obtain the "generic" spectral behaviour. Let us examine the 
Ihs of the last relation. To this aim, we denote 

K c :=4(|t n | 2 + |t 22 | 2 + |t 12 | 2 + |t 21 | 2 ) 
fr s :=4(l + |t n t 22 -t 12 t 21 | 2 ) 
L c (6 1} 9 2 ) : = 8 [-3? ((£11*21 + t 12 t^)e m ) + 3? ((i 22 *2i + hit 12 )e l92 )} 
L(6 U 9 2 ) := 83? {tnt^e^-^) + 83? (t 12 & 01 ^) 

which allows us to write the coefficient V 2 = V 2 (ak, 6>i, 9 2 ) as follows 

V 2 (x,6 l} 6 2 ) = -K c cos 2 x + K s sin 2 x + cosx ■ L c (6 1 ,6 2 ) + L(9 1 ,9 2 ) ■ 

To examine the spectral asymptotics, the following functions, 

V 2 + (x) :=max{V 2 (x,e 1 ,e 2 )\e 1 ,e2 e (-vr,7r]} , 
V 2 -(x) :=mm{V 2 (x,e u 6 2 ) \ 6 U 6 2 G (-tt.tt]} , 

will be essential, since the spectral condition can be expressed, up to an error of order 
of O (k" 1 ), by the inequalities, 

V 2 + (ak)>0 A V 2 '(ak)<0. (16) 
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As we will see below, the following two constants will play an important role: 

L+ := 8 max {ft ({tut^ + t 12 t^)e Wl ) + ft ((t 22 t^ + t^t 12 )e i62 ) + 
+SR(t u h~ 2 e i{ei - e2) ) +^(t 12 t 2 - 1 e i ^ +d ^) \ 9 1 ,9 2 G , 
L~Z :=8min{ft(t n &( ei -^) + ft (t^e^ 1 ^) | 9 U 9 2 G (-7r,7r]} ; 

it is easy to see that = —8\tnt 22 \ — 8 1 ^12^21 1 - 

2 

With this preliminary we are going to formulate and prove a claim which will be 
useful not only here, but also at other places further on. 

Proposition 6.1. Let 

F(9 1 , 2 ) = ft (A^ 01 ) + ft (A 2 e W2 ) + ft (Age^ 1 ^) + ft (A^ 1 ^ 2 )) , 

where the coefficients Aj, j = 1,2,3,4, are independent of8i,8 2 . Then the range of this 
expression, F := {F(6i,6 2 ) | 6\,6 2 G (—ti,tt]} ; is an interval which is non- degenerate if 
and only if there is an index j G {1, 2, 3, 4} such that Aj ^ 0. 

Proof. Since F is continuous and (— 7r,7r] 2 is connected, F is an interval. To finish the 
proof it remains to check that a constant C G K such that 

F(9 U 9 2 ) = C for all (9 U 9 2 ) G (-vr, vr] 2 (17) 

exists if and only if Aj = for all j = 1, 2, 3, 4. 

Consider a fixed G K and a number A G C such that argA = a, then 
ft(,4e i9 ) = ft(|A|e ia e ie ) = |A| cos(a + 9). We apply this idea to rewrite ([T?D as 

I Ai| cos(cti + 6 l i) + |v4 2 | cos(a 2 + 2 ) + |A 3 | cos(a 3 + 9 1 + 9 2 + \A 4 \ cos(a i + 9 1 -9 2 ) -C = 

for all (9i,9 2 ) G (— 7r,7r] 2 . It is easy to see that the 5-tuple 

1, cos(«i + #1), cos(a 2 + 9 2 ), cos(a 3 + 9 1 + 2 ), cos(a 4 + #i - # 2 ) 

is a linearly independent set of functions on [0, 27i~) 2 , therefore (fTTj) is satisfied if and 
only if A 1 = A 2 = A 3 = A± = C = 0. □ 

To solve the spectral conditions ( JT6|) . we need to know basic characteristics of the 
functions V 2 + (x) and V 2 ~{x) involved in it. 

Proposition 6.2. The functions V 2 (x) and V 2 (x) have the following properties: 
(i) Both V 2 + (x), V 2 (x) are it -periodic and satisfy 

V 2 + (^-x)=V+(x) 1 V 2 -(^--x)=V 2 {x). 
(ii) Function V 2 ~(x) is increasing in [0, f] and 

V 2 -(0)<0, V 2 -(|) =K S + L~„. 
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(Hi) It holds V^ + (0) = — K c + Lq , and there is a number Xq G [0, |] such that 

— V 2 is increasing in [0,xq], 

- V 2 + (x) > for all x G [x , |] . 

(iv) If at least two entries ofT are nonzero, then V 2 + (x) > V 2 (x) for all x G (0, f). 

Proof of this proposition is technical and slightly long, and can be found in the 
appendix. 

Using the above result we can characterize the spectrum. 
Proposition 6.3. If at least two entries ofT are nonzero, then the following holds true: 

• If Lq > K c and L^_ < —K s , there is a ko > such that the interval [k^, +oo) is in 
the spectrum. 

• If Lq < K c , the spectrum has asymptotically gaps of the form 

fn 2 ir 2 2bnir , n 2 ir 2 2bnir M , A 

— + ° 1 > — + — + 01 > 

\ a 1 a 1 or a z J 

where b G (0, f) is the number uniquely determined by the condition V 2 + (b) = 0. 

• If L^ > —K a , the spectrum has asymptotically gaps of the form 

2 

where c G (0, |J is uniquely determined by the condition V 2 (| — c) = 0. 

Proof. The argument is based on Proposition 16.21 (i)-(iv). First we notice that if 
Lq > K c , then V 2 + (x) > for all x G [0, |] according to (^zzj, and thus for all x G [0, +oo) 
due to (i). Similarly, if L^_ < —K s , then V 2 (x) < for all x G [0, |] according to fmj, 
and thus for all a; G [0, +oo) due to (i). Consequently, for both Lt > K c and L z < —K s , 
the asymptotic spectral condition (|T6|) is satisfied for all k. 

If Lq < K c , then by virtue of (ii) there is a 6 G (0, |) such that V 2 (x) < on [0, b) 
and V 2 (x) > on (b, |). Then the first inequality of (TTdT) is not satisfied on (0, and 
we infer from (i) that it is not satisfied on each interval (— — -, — + -). 

In the same vein, the inequality > — iT s in combination with (Hi) implies 

2 

existence of a c G (0, |) such that V 2 (x) < on [0, | — c) and V^~(x) > on (| — c, |) . 
The second inequality of (fTBT) is then not satisfied on — ^ , ^ + £ ), and by (i) it is 
not satisfied on each interval ((n + |) - — - , (n + |) - + - ) . 

Finally, it follows from ^roj that if the matrix T has at least two non- vanishing 
entries, we have b < | — c and the gaps cannot overlap asymptotically. □ 
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6.2. Further particular cases 

Proposition 16.31 applies if at least two entries of T are nonzero. If this is not the case, 
the spectral condition simplifies significantly and one of the following situations occurs: 

• S12 = 0. With such boundary conditions, the lattice is decomposed into separated 
edges or pairs of edges, and consequently the spectrum is pure point. As above in 
similar cases it is infinitely degenerate, of course, but the contribution from each 
edge pair has the usual semiclassical behaviour: the number of eigenvalues not 
exceeding k 2 has Weyl asymptotics with the leading term -. 

2 2 

• s±2 7^ and T = 0. The spectrum contains points ZL -§ L and bands of asymptotically 

2 2 2 2 

constant width in the vicinity of the points may or may not lie in the 

bands. The gaps grows linearly as n — > oo. 

• S12 7^ and T / 0. The spectrum asymptotically consists of bands located in 
the vicinity of the points + ~ arccos jrpjfj^ J j where t is the nonzero entry of T. 
They are of asymptotically constant width while the gaps are linearly growing as 
the band index n — > oo. 

It remains to deal with the special cases Lt = K c , = —K s left out in Proposition ^. 31 

2 

We will not go into much detail here, but it is worth to solve this situation in the 
particular case of scale-invariant coupling, i.e. for S = 0. It turns out that: 

• If S = and Lq = K c , then there are no gaps around 

• If S = and = —K s , then there are no gaps around (n + |) 2 

7. The case of m = 3 

In this case we have a Hermitean S and complex numbers tj determining 

Sl3 \ 

■S33 

-h 1 / 

In contrast to the previous case, there is no problem with the renumbering invariance, 
because the edge numbering can be changed if necessary by a trivial rotation of the 
lattice. A direct calculation of the determinant in (fTUj) yields the spectral condition, 



B 



( 


1 








tl 


\ 







1 





tl 












1 






\ 














/ 



A = - 







S12 




S12 








$23 


\ 


-h 


-h 



V 3 ■ k 3 + V 2 ■ k 2 + V x ■ k + V = , 
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where Vj, j = 0, 1, 2, 3, are expressions depending on ak, on the entries of the matrices 
S, T, and on the quasimomentum components 9i,9 2 , given explicitly by 

V 3 = 4sinak [(l + |ti| 2 + \t 2 \ 2 + \t 3 \ 2 ) cosak + 23?(t 1 t 2 ~e iei ) - 23?(t 3 e ie2 )] , 

V 2 = 4 cos 2 ak [-(sn + s 22 ) (1 + |t 3 | 2 ) - s 33 (\t 1 \ 2 + |t 2 | 2 ) + 23? ((s 13 ^ + s 23 t^)t 3 )] 
+ Asm 2 ak [s u \t 2 \ 2 + s 22 |^i| 2 + S33 - 23? («i 2 M 2 )] 
+ 8cosak [3? ((-S12 + S23M3 + S13M3 - si 2 |t 3 | 2 - s 33 tih)e iei ) 

+3? + s 22 )t 3 - s^h - s£t 2 ) e id2 )] 

+ 83? ((s 12 t 3 - Wi)^ 01 ^) + 83? ((s 12 T 3 - s 13 r 2 )e^-^) , 

V 1 = Asinakcosak [(\s l2 \ 2 - s n s 22 )(l + |t 3 | 2 ) + (\s 13 \ 2 - Sns 33 )(l + |t 2 | 2 ) 

+ (k 23 | 2 - S 22 S 3 3)(1 + l^l| 2 ) + 

+23? ((si 2 s 33 - s^s^hh) + 23? ((S23S11 - s 13 s^)ht 3 ) + 23? ((si 3 s 22 - si 2 S23)*i*3)] 
+ 8 sin ak ■ 3? ((S13S23 - s 33 s 12 )e l61 ) 

+ 8sinak ■ 3? [((si 2 s 23 - s 22 sli)ti + (s 12^ ~ s u s^)t 2 + (s u s 22 - |si 2 | 2 )t 3 ) e 102 ] , 
V = -4 sin 2 ak ■ det S . 

The case m = 3 is probably the most interesting one. While in the one-dimensional 
case we have bands and gaps which are asymptotically either of constant width or 
linearly growing with the band index, a square lattice with m = 3 exhibits a considerably 
richer spectral behaviour. The band are here of two types, mutually interlaced: 

2 2 

• "even" bands in the vicinity of IL -j- with n G N, 

• "odd" bands in the vicinity of (n + |) 2 ^ with n G N, 

In this section we focus on the generic situation, \ti\ ^ \t 2 \ or \t 3 \ 7^ 1, when we 
will be able to describe several interesting asymptotic types. The other possibility, 
\ti I = \t 2 \ A \t 3 1 = 1, will be briefly commented at the end. 

Let us start with the spectral condition in the form V 3 = 0(k^ 1 ), i.e. 

4sinaA; [(l + |tx| 2 + |t 2 | 2 + |t 3 | 2 ) cosaA; + 23?(t 1 t 2 V f?1 ) - 23?(t 3 e ie2 )] = G(k~ 1 ) . (18) 

The Ihs of the last equation should be close to zero for large k, which can be achieved 
either for small absolute value of sin ak or for small absolute value of the expression in 
the brackets. These possibilities refer to the "even" and "odd" bands mentioned above, 
respectively; they will be described in detail below. Unless stated otherwise, we suppose 
that |fx| ^ \t 2 \ V |i 3 | ^ 1. 
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2 2 

7.1. Generic case: asymptotically constant spectral bands around 

We start with the bands corresponding to small values of | sinafc|. Let us consider the 
spectral condition in the form V 3 + ^ = 0(k~ 2 ), divide it by four and rewrite it as 

sinaA; [(l + |*i| 2 + \t 2 \ 2 + \t 3 \ 2 ) cosak + 2$l(t 1 t 2 ~e Wl ) - 23?(t 3 e i92 ) 
sin ak 



k 

cos 2 ak 



(s n \t 2 \ 2 + s 22 \t 1 \ 2 + s 33 - 23? (si 2 tit 2 )) 



[-(sn + s 22 ) (1 + |t 3 | 2 ) - s 33 (|^i| 2 + M 2 ) + 23? ((si 3 ti + s 23 t 2 )t 3 )] 

- [3? ((-su + s^hh + suhta - s 12 \t 3 \ 2 - s^hh)^ 01 ) 

+3? (((s u + s 22 )t 3 - slit! - s£t 2 ) e ie2 )] 

- ((auta - ^i)e i(9l+e2) ) - ^3? (( fll2 £ - s l3 T 2 )e^~^) + 0(fc- 2 ) . (19) 

We restrict our considerations to the values of k such that |cosa&| > c where c is a 
constant satisfying ^q^pqrjg^^ < c < 1; recall that such c exists in view of the initial 
assumption \t\\ ^ \t 2 \ V \t 3 \ ^ 1. Due to this restriction, the absolute value of the term in 
the brackets at the Ihs of (Tl9|) can be asymptotically estimated from below by a positive 
constant, whence it is easy to see that high in the spectrum we have sin ak = 0(k~ 1 ). 
We set 

nir d 
k = — + - , 

a n 

then sinafc = (-l) n ■ f + C(n" 3 ), cosaA; = (-l) n + 0(n- 2 ) and § = ^ + 0(n- 2 ). 
We substitute from these relations into (Tl9|) . divide both sides by the expression in the 
brackets at the Ihs and after elementary manipulations we arrive at 

— • [(1 + |ti| 2 + \t 2 \ 2 + |t 3 | 2 ) + 2(-l) n $t(t i r 2 e^) - 2(-l)"3?(t 3 e ie2 )] = 



n 

a 



— [-(sn + s 22 ) (1 + |t 3 | 2 ) - s 33 (|ti| 2 + |t 2 | 2 ) + 23? ((s 13 ti + s 23 t 2 )t 3 )] 
l) n a 



^—^T- I K ((~ s i2 + Wztxh + s 13 t 2 t 3 - s 12 \t 3 \ 2 - s 23 M 2 )e i(?1 ) 
+3? (((s u + s 22 )t 3 - s^h - s^t 2 ) e W2 )} 



2a „ „ . „ fl lfl ^ 2a 



- — « ((a 12 t 3 - ^i)e i(ei+e2) ) 3? (( Sl2 ^ - s 13 r 2 )e^^) + 0(n~ 2 ) . 

We observe that all the terms ( — l) n can be "absorbed" into d\ and 9 2 by the shift 
(6i,9 2 ) i — ^ (^i + nir,6 2 + mr), and thus neglected. Therefore the asymptotic condition 
can be written in the form 

1 t 

for 

W 3 (T, 1; 2 ) := (1 + |t!| 2 + |t 2 | 2 + |t 3 | 2 ) + 2(-l)"3?(t^e iei ) - 2(-l)"3?(t 3 e ie2 ) , 
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w 2 (s,T,e 1 ,e 2 ) : = -( Sll + s 22 ) (1 + |t 3 | 2 ) - s 33 (|til 2 + |t 2 | 2 ) + 23? ((s 13 r 1 + s 23 r 2 )t 3 ) 

+ 23? ((-s 12 + s^titg" + sis^ts - si 2 |t 3 | 2 - s 33 t!^)e iei ) 
+ 23? + s 22 )t 3 - slit! - si£t 2 ) e i02 ) 

+ 23? ((s 12 t 3 - WzhW^ 9 ^) + 23? ((s 12 T 3 - s 13 F 2 )e^-^) . 



Now we denote 

-w 2 (s,T,e 1 ,e 2 ) 



01,02 e (-7T,7r] 



W 3 (T,0 h 2 ) 

note that T> is an interval which is bounded due to the condition \t\\ ^ \t 2 \ V |t 3 | ^ 1. 
The spectral bands are then described in terms of T> as follows, 

11 7T d 

k = — + — + 0(n' 2 ) , deV, (20) 
a nix 

which yields the allowed energy values 

fc 2 = ^ + 2- + 0(^), dev. 

a 1 a 

Consequently, these bands are of asymptotically constant width as n — > oo. 
Furthermore, it is easy to show that if 

(a) T = and Si 2 = 0, or 

(b) t\t 2 = 0, t 3 = and 5* is diagonal, 

then the fraction ^^(r'gfg^ 2 ' 1 is independent of (6*i, 9 2 ), and thus the interval V shrinks 
to a point. These cases need a special approach. 

One can check that (b) refers to the situation when the lattice is decomposed either 
to individual edges (if T = 0) or to pairs of edges (if T ^ 0), with infinitely degenerate 
eigenvalues in the vicinity of (—J . The case (a) is much more interesting and we will 
analyze it below. 

7.2. Spectral bands around shrinking as n~ 2 

Under the assumptions of the case (a) above the spectral condition simplifies to 

4/c 3 sin ak cos ak + Ak 2 \— (sn + s 22 ) cos 2 ak + s 33 sin 2 ak\ 

+ Aksinak [(|si 3 | 2 + |s 23 | 2 - s n s 22 - (s n + s 22 )s 33 ) cosak + 23? (siss^e 101 )] 

— 4 sin 2 ak ■ det S = . 



The spectrum is given by fl20l) with T> = (s n + s 22 )}. Since the latter is a one-point 
set, the spectral band is determined by the next term in the expansion. We set 
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and performing a calculation similar to that of Sec. 17.11 we arrive at the expression for 
d! . Surprisingly, dl itself behaves like n _1 , i.e. the last term above is of order of 0(n~ 3 ), 



d! a(s u + s 22 ) 



n 2 n 3 n 3 



( s 2 + s 2 \ 
[sn + s 22 ) + al 11 2 22 +2s llS22 + \s 13 \ 2 + \s 23 \ 2 \ - 2a$l (s 13 s^e l61 ) 



Thus k belongs asymptotically to the spectrum if 



_ n-n t gn + s 2 2 , ajsn + s 22 ) 
a nir n^ir^ 



-{sn + s 22 ) + a I " u * 522 + 2sns 22 + |si 3 | 2 + |s 23 | 



2a 2 

+ ^rl s n + s 2 2 | • IS13S23I cos^ + 0{n 4 ) for ■& e [0, 2vr) . 



n 3 7r 3 

The left and right endpoint of the spectral band are then given as k\ and k 2 R with = tt 
and $r = 0, respectively, hence the band width behaves like 

k 2 R- k l = i k R ~ k L ){k R + k L ) = 8-^-=-|sn + s 22 | ■ |si 3 s 23 | + 0{n~ 3 ) 

in the asymptotic regime, n — > 00. 

Note that there are two exceptional cases here, namely: 

• S13 = or s 23 = 0. The spectral condition is independent of 81, 9 2 , thus it has only 
point solutions. 

• sn + s 22 = 0. The spectral condition factorizes, sin ak [cos ak + O mm] = 0, 
producing no (true) band in the vicinity of 



7.3. Generic case: linearly growing spectral bands around (n + |) 2 ^ 



Now we proceed to the "odd" bands, i.e. those corresponding to asymptotically small 
values of the expression |(1 + |ti| 2 + |t 2 | 2 + |t 3 | 2 )cosa£; + 2$l(t 1 t 2 ~e Wl ) - 2^(t 3 e i02 )\. 
Their structure can be derived directly from the spectral condition (fl8|) . we just replace 
for the sake of simplicity the term 2$l(t 1 t 2 ~e i91 ) - 2$t(t 3 e id2 ) by 2(|M 2 | + |t 3 |)cosi? in 



analogy with what we did in Section 15. 2\ cf. Observation 15.11 

4sinafc [(l + \h\ 2 + \t 2 \ 2 + |t 3 | 2 ) cosafc + 2(|tit 2 | + \h\) cost?] = O(k^) 

Suppose that |tit 2 | + |t 3 | 7^ (the case t\t 2 = t 3 = is postponed to Sec. 17.41 below), 
then we may divide the equation by 8(|£]t 2 | + |t 3 |) obtaining 



sin ak 



1 + 


tl 


2 + 




2 + \ts\ 2 




2(1 


tit 2 \ 


+ 


l*8|) 



cos ak + cos $ 



0(k~ L ). (21) 



Since the situation when | sina/c| is very small has been treated in the part devoted to 
"even" bands, here we adopt the premise that | sina/cj is large enough. More precisely, 
we restrict our considerations to such values k that | sinafc| > c, where c is a constant 
satisfying < c < 1, the exact value plays no role. Then (l2~Tj) can be rewritten as 



1 + 


tl 


2 + 


1*2 


2 + |* 3 | 2 


2(|tit 2 | 


+ l*s|) 



cos "d = j : : — p- cos ak + 0{k 



-i> 
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Such a condition has a solution for some parameter $ if the modulus of the rhs rhs is 
smaller than one, i.e. for 



cos ak\ < - u i2 u I u , 2 + • (22) 







1 + 


*1 


2 + l*2l 2 + 


ts\ 2 



In analogy with Sec. 15.21 we infer that the spectral bands are neigbourhoods of the 
points [(n + |) ^] 2 . We set k = (n+ |) - + c?, thus | cosafc| = | sina<i| and \ of the 
n-th band solution is 0{n~ l ). Then ( 1221) can be rewritten as 



smad\ < - + 0(n x ) . 



2(|M 2 | 


+ 


N) 


1 + 


*i 


2 + 


1*2 1 


2 + 


h 


2 



Recall that it follows from the initial assumption \b%\ ^ | *2 1 V | *3 1 7^ 1 that the first 
term at the rhs is smaller than one, cf. Observation 15.11 hence it makes sense to set 



A := arcsin p^p^^j^j^p £ (0, |). The spectral bands are then given by 

1 \ 2 7T 2 nil , „„ , . / 1 \ 2 7T 2 727T . „ . . 

in the high-energy asymptotics, n — > 00. 



7.^. Asymptotically constant spectral bands at (n + ^) ^ 

Let the matrix T satisfy tyt 2 = i 3 = 0. We may suppose without loss of generality that 

V2 

k 



t 2 = 0. The substitution t 3 = t 2 = into the spectral condition V3 + = C(& 2 ) gives 



/ 1 .ox cos ci/u / 1 ,o\ sin dk / 1 1 9 \ 
smafccosafc (1 + |*i j j = — - — [s u + s 22 + s 33 |*i| j — (s 22 |*i| + s 33 ) 

+ [X (s 12 e^) + K (s^e?*)] + ?H (f^e 1 ^) + 0(AT 2 ) . (23) 

We proceed in a way analogous to Sec. 17.11 setting k = [n + |) - + 5 and substituting 
the expansions 

cosaA; = -(-l) n ■ aS + O (S 3 ) , sinafc = (-1)" + O (S 2 ) , i = — + e>(n~ 2 ) 

into (I23p . which yields after a simple manipulation an expression for S, 
1 



5 = 
nil 



S22 




2 + s 33 


1 + 


1*1 


2 



+ TTir^l^( t « i(ei+e2) ) 

1 + \ti r 



0(n 



-2\ 



Hence the spectral bands are determined by the values of satisfying 

, / . 1\ 7T 1 S 2 2|tl| 2 + S 33 . 1 2|t 1 S 23 | q . _2n 

fc= n + - • : — ^ • ; — [77 cos v + Cm 

V 2/ a nir 1 + 2 nvr 1 + U 2 V ; 
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so that 



1 \ 7T 



\ 2 / a z a 1 + a 1 + 

giving rise to asymptotically constant-width bands. 

Note that the leading "band-producing" term may collapse to a single point which 
happens if t\ — or S23 = 0. The first named situation will be discussed in detail in 
Sections 17. 5H7. 71 the second one will be omitted, since the analysis is in principle similar. 

7.5. Spectral bands at (n + |) 2 ^ shrinking as n _1 

The bands in the vicinity of (n + ^ may shrink as n' 1 provided the matrices 5* and 
T satisfy certain conditions in addition to t 2 = t 3 = 0, namely t\ — and S13S23 7^ 0. 
We will demonstrate this fact in the first case, £1 = 0. 

We proceed as in Sec. I7.4[ but we take one more term in the expansion of k, i.e. 

we set 

k= ( n+ l)i + «i + s, 

\ 2/ a mi 

and substitute into the spectral condition in the form V3 + ^ + p- = 0(k~ 3 ). A 
calculation leads to the expression 5 = \ [— ^ + ^§9R ((si3S23)e l9l )j in the leading 
order, hence 

1^ 7T i S 33 S33 , 2a ,_ _ , , ,„ f _ 3 



fc= n +o - + — - TT^" + ^^l s i3S 23 | costf + 0(n _d ) , ■& 6 [0, 2tt) . 

The band edges correspond to = 7r, 0, respectively, and the band width equals 

g 

— \s 13 s 23 \ + 0(n~ 2 ) , 
mi 

i.e. the bands shrink asymptotically as n~ l unless s 13 = or s 2 3 = - these cases will 
be discussed in the following sections. 



7.6. Spectral bands at (n + |) 2 ^ shrinking as n 
If T = and S23 = 0, the following spectral band asymptotics can be derived, 

1 \ 7T , S33 S33 i C a S33 i „ , q , /ri/ __4\ 



-2 



k =\ n +n)- + — ~T^- + -3--3- HT^I C ° S ^ + °( U ) > * G I ' 2?f ) 

where C depends only on S and a. This determines bands of the widths 

4 

l s 33 s 12 \ + (D(n~ 3 ), 



n 2 ir 



i.e. they shrink asymptotically asn 2 providing that S33 7^ and s\ 2 7^ 0. 

Remark 7.1. If T = 0, S23 = and s\ 2 — 0, the spectral condition is independent of 
#1, 6*2, thus the spectrum is pure point. 
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7.7. Spectral bands at (n + |) 2 ^ shrinking 



as n 



Let T = and s 23 = s 33 = 0. Then the spectral bands are characterized by the following 
asymptotical condition, 



7T a 2 |si 3 | 2 s 22 
n + — 



3|si 3 | S22 o_ 

A 



(I 



n° 



7T° 



n 1 



2tt 3 



rr 



2|-Sl3| 2g 22 

2k a 



S12I COS-i? + 



with d e [0, 27r), thus the band widths are given by 



a 2 |si 3 | 2 s 2 2 



S12I cos-# + 0{n 



n 3 2tt 3 

i.e. they shrink asymptotically as n~ 3 providing s 12 7^ 0, S13 7^ and s 22 7^ 0. 
Remark 7.2. Let T = 0, s 23 = and s 33 = 0. Then: 

• If S\2 = 0, the spectral condition is independent of the quasimomentum components 
0i, 02 so the spectrum is pure point. 

• If S13 = or s 22 = 0, the spectral condition can be factorized in the form 
cos ak ■ [sinak + O (-)] , and consequently, the band in the vicinity of (n+ |) ^ 
collapses into a single point. 

7.8. The particular case \t\\ = |t 2 | A |t 3 | = 1 

This situation will not be treated in detail, but we will comment on the prominent 
case 5* = which corresponds to the scale-invariant vertex coupling. If 5* vanishes, the 
spectral condition simplifies to V 3 = which can be rewritten as 

4sina£; [(2 + 2|ti| 2 ) cosaA; - 2|ti| 2 cos$i - 2cos$ 2 ] = , 

where and ?9 2 are properly shifted 0\,02- The choice $1 = $ 2 = ak sets the expression 
in the brackets to zero, therefore the spectrum contains the positive halfline. 

8. The case of m = 4 



Finally, we pass to the case which is generic from the viewpoint of the boundary 
condition ([1]) when the matrix B is regular. Following the discussion in the opening 

we set 

/ Sll «12 -sra \ 

51 2 S22 S23 S 2 4 

513 ^3 S 33 S 34 
\ Sl4 S 2 4 S 3 4 S44 J 

There is obviously no problem with the renumbering the lattice edges. A direct 
calculation of the determinant in (fTUj) leads to the spectral condition 



B 



( 1 \ 

10 

10 

1 



.4 



V 4 ■ k 4 + V3 • k 3 + V 2 ■ k 2 + V 1 ■ k + V = 
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where Vj, j = 0, 1,2,3,4, are expressions depending on ak, on entries of the matrix S, 
and on the quasimomentum components 61,62, given be the formulae 

V4 = —4 sin 2 ak , 

V 3 = 4 sin ak [(% + s 22 + s 33 + s 44 ) cos ak + 23R (si 2 e iei ) + 23ft (s u e W2 )] , 

V 2 = 4 cos 2 ak (\s 13 \ 2 - sns 33 + |s M | 2 - s n su + \s 23 \ 2 - s 22 s 33 + |s 24 | 2 - s 22 s A *) + 
+ 4sin 2 a/c (s u s 22 - \s 12 \ 2 + s 33 s 44 - \su\ 2 ) 
+ 8cosak [-(333 + s 44 ) ■ 3ft (si 2 e i91 ) + 3ft ((s 13 s^+ si 4 s5l)e Wl ) 
-(six + S22) • 9? (W 2 ) + « ((slis 14 +» 24 )e i02 )] 
+ 83ft [( Sl4 s^ - Sl2 s 3 4)e i(ei+92) ] + 83ft [(s 13 sii - s^e 1 ^- 02 )] , 

Vo = —4 sin 2 afc • det S . 

The remaining term V\ will not be needed - it is sufficient to know that it is bounded 
with respect to the parameters. The spectral condition can be written as 

sin 2 ak . [( Sn + S22 + S33 + S44 ) cos a k + 23ft (si 2 e iei ) + 23ft (s 34 e i92 )] 

H p — (I s i3| 2 - sns 33 + |si 4 | 2 - sns 44 + |s 23 | 2 - s 22 s 33 + |s 24 | 2 - s 22 s 44 ) 

+ [-(ass + ^44) • H (, 12 e iei ) + K ((s 13 ^ + Sl4 sil)e i91 ) 

-(su + s 22 ) • 3? ( S34 e ie2 ) + 3ft ((5Hs 14 + <^ 24 )e ie2 )] 

+ ^ [((*i4*5S - s 12 s 34 )e i(eiHN?2) ) + 83? ((si355I - s^e 1 ^-* 3 ))] 

sin ci/u / 1 1 2 11 2\ 

H ^ — {sns 22 -\si 2 \ + s 33 s 44 - |s 34 | J 



k 2 

O (at 3 ) . 



(24) 



Since the r/is is of the order of 0(k~ 2 ), it follows that the values of k that solve the 
spectral condition are asymptotically close to the points n G N. With this fact in 
mind it is convenient to express k introducing d such that 

k = 1 . (25) 

a nit 



Substituting this into (1241) . we get after a simple manipulation, 
d 2 - d ■ [s n + s 22 + s 33 + s 44 + 2 ■ (-1)" ■ 3ft (s 12 e^) + 2 ■ (-1)™ • 3ft (s 34 e ie2 )] 

(\ |2 1 I |2 1 I |2 I I |2 \ 

— U s 13| — -h — SnS 44 + |S 23 | — S 22 S 33 -+- |S 24 | — S 22 S 44 J 

- 2 ■ (-1)™ • [-( S33 + s 44 ) • 3ft (s 12 e idl ) + 3ft ((S13SE+ sus^e* 1 ) 
-(sn + s 22 ) • 3ft (s 34 e ie2 ) + 3ft ((sii-sw + ^ 24 )e i92 )] 
- 23ft [(s 14 ^i - si 2 s 3 4)e i(ei+(?2) ] - 23ft [(s l3 s^ - s 12 sE)e^-^] = O {n' 1 ) . (26) 
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We observe that the last relation has the structure 

d 2 + 2p ■ d + q = O (n- 1 ) (27) 

with the quantities p, q depending only on S, 6\, 62 and the sign of n. A value k 2 > 
with k of the form (|25|) belongs to the spectrum if there are 61, 62 such that the above 
condition is satisfied. The solutions of (1271) are given by 



d lfl = -p± ^p 2 -q + O (n- 1 ) . (28) 
Proposition 8.1. It holds max{p 2 — q \ 6\, 62 £ (— 7T, n]} > 0. 
Proof. Taking the value of q from (126|) and p 2 arranged into the form 

P 2 = \ ( S H + S 22 - S33 - S44) 2 + (S11S33 + S11S44 + S22S33 + S22S44) + 

(-1)™ ■ + S22 + s 33 + s M ) [n (s 12 e^) + to (s 34 e i9 >)] + [K (^i 2 e iei ) + 3* (s 34 e ie2 )] 2 , 
we find 

P 2 ~ q = \ («n + s 22 - s 33 - ^44) 2 + [3* (si 2 e Wl ) - 3* (s 34 e ie2 )] 2 

+ (-l) n • 3* [((s u + s 22 - s 33 - s 44 ) s X2 + 2 (s 13 si£ + s 14 sii)) e iei ] 
+ (-1)" • 3? [((-an - s 22 + s 33 + s 44 ) s 34 + 2 (s7Js 14 + S2is 24 )) e i9a ] 
+ |si 4 | 2 + h 3 | 2 + 29* (s u sr^ m+e2) ) + \si 3 \ 2 + \s 2 ,\ 2 + 23* (auSSe 1 ^-^) - (29) 

It obviously holds 

|s 14 | 2 + hsl 2 + 23* {s^e^ 9 ^) > (\s u \ - \s23\f > , 
|s 13 | 2 + \S24\ 2 + 23* (si 3 » i(9l - e2) ) > (|si 3 | - |s 24 |) 2 > , 

and 61,62 can be chosen so that the terms on the second and third line of (129]) are 
non-negative. Therefore 

maxjjo 2 — q \ 61,62 & (— it, n]} 

> 7 ( s n + s 2 2 - s 33 - s 44 ) 2 + (|s 14 | - |s 23 |) 2 + (|si 3 | - |s 24 |) 2 > , 

which we have set up to prove. □ 

Corollary 8.2. If max-jjo 2 — q\0i,02 £ (—11,71]} = 0, then S = diag(sn, s 22 , S33, S44) 
and s u + s 22 - s 33 - s 44 = 0. 

Proof. It follows from the last part of the previous proof that the premise requires 
Sn + S22 — s 33 — S44 = 0. If we substitute this into rhs of ( )29|) and make the same 
estimate as in the proof, we infer s 13^3 + s 14S24 = 0, S13S14 + S23S24 = 0. Then 6\, 62 can 
be obviously chosen such that p 2 — q = (\sn\ + | S23 1 ) 2 + (| Sx3 1 + l s 24|) 2 , hence necessarily 
s u = S23 = S13 = s 24 = 0. Now d29]) gives p 2 - q = [3* (si 2 e iei ) - 3* (s 34 e i92 )] 2 which 
does not exceed zero only if s\ 2 = S34 = 0. □ 
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8.1. Point spectrum for a diagonal matrix S 

If m = 4 and all the off-diagonal entries of S vanish, then the boundary conditions 
(JIJ) describe the system of decoupled edges of the length a. Such system has a point 
spectrum which is not difficult to find. We have two edges types: 

• "horizontal", with the boundary conditions i/)'(0) = 32^(0), V ; '( a ) = ~~ s n' l P{ a )y 

• "vertical", with the boundary conditions ip'(0) = 544^(0), ip'(a) = —s 33 ip(a). 

The corresponding spectral conditions are 

[Sl+2j,l+2j + kS2+2j,2+2j) COS ka = (k — Si +2 j t \ + 2 jS 2+2 j, 2+2 j) SHI &d , j = 0, 1 . 

Their solutions are thus 

k = mr + aictga^,^ + s 1+2i , 1+2i + Q , 2) 
a a nn 

giving rise the eigenvalues, or flat bands of the lattice Hamiltonian, 

k 2 = ( nn ) 2 | 2 w7r ajctgg2+2j,2+2j | / arctgg 2+2j - 2+2 A 2 | g s^i+gj | Q / - 2 \ 
\ a J a 2 \ a J a 

for j = 0, 1 corresponding the two edge orientations. 

8.2. A general matrix S 

Let us next consider the following sets related to solutions to condition (128]) . 

V lt2 = {- P ± vV-<? I Ox, 02 e (-7i, tt]} n R , 

It follows from Proposition 18. II that T>\ ^ 7^ P 2 . We will show now that if one of them 
contains more than one point, then both contain a non-degenerated interval. 

Proposition 8.3. Let one of the expressions —p+^/p 2 — q, —p—yjp 2 — q be independent 
of 62), then the same is true for the other one. 



Proof. Let us suppose that — p + ^p 2 — q = c with c £ R holds for (#1, 2 ). In such a 
case \/p 2 — q = p + c, hence — q — 2pc = c 2 , in other words, the expression — q — 2pc is 
independent of (#1,6^). Since 

— q — 2pc = ( I S13 1 2 — S11S33 + I S14 1 2 — S11S44 + I S23 1 2 — s 22 s 33 + \s2i\ 2 ~ S 22 .S44) 

+ 2 • (-1)" • [-( S33 + s 44 ) • 3? (s l2 e^) + K ((s 13 s£ + Si^K* 1 ) 

+ s 22 ) • 3? ( S34 e ie2 ) + SR ((slis 14 + ^ 24 )e ie2 )] 
+ 29* [(5 14 S2i - s l2 s 34 )e [ ^ + ^] + 2U [(s 13 sE - s l2 s^)e^^] 
+ c ■ [su + S22 + s 33 + s 44 + 2 • (-1)" • H (s 12 e iei ) + 2 ■ (-1)" ■ 3? {s 34 e i6 *)} , 
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we easily find that also the expression 

(-l) n • 3? [((-(S33 + S44) + c) • si2 + s 13 s^ + s u s^) e Wl ] 

(-l) n ■ 3? [((-(sn + s 22 ) + c) ■ s u + s^s u + <^s 24 ) e i9a ] 

+ & [(si 4 s5i - s 12 s 3 4)e i(ei+e2) ] + 3ft [(siaSM - sia^e^ 1 -^)] 

should be independent of (6\, 62) ■ It follows from Proposition 16 . 1 1 that this is true if and 
only if the whole expression identically equals zero. This in turn means that the Ihs of 
the condition (127|) is independent of 61,82, and consequently, both roots of the equation 
d 2 + 2pd + q = are independent of 9\, 6 2 . □ 

Proposition 18.31 in fact says that only two situations are possible, either both sets 
T>\, T>2 are non-degenerate intervals, or each of them contains a single element only. 

Proposition 8.4. The sets T>\, T>2 are single- element sets iff at most one pair of off- 
diagonal elements of S is nonzero, and moreover, s 12 = S34 = 0. 

Proof. T>i, T>2 are one-element iff both — p + y / p 2 — q, —p — y / p 2 — q are independent 
of 61,62, and this obviously holds iff both p, q are independent of 61, 6 2 . With regard 
to the definition of p, q, we have 

S12 = , s 34 = , S13S27 + S14S2I = , slisi4 + S2IS24 = , 
•S13S2I = , S14S2I = . 

Hence s±2 = 0, S34 = 0, and it is easy to see that at least three elements from the set 
{S13, S14, S23, s 24 } have to vanish as well. □ 

Corollary 8.5. The sets T>\, P 2 are one-element sets iff the lattice decouples either into 
separate edges of the length a, or into identical copies of L-shaped pairs of edges. 

Proof. See the previous proposition. If all the off-diagonal elements of S vanish, the 
system decouples into separate edges and we return the the situation considered above. 
The second possibility is that just one of the numbers S13, S14, s 2 3, s 2 4 is nonzero, then 
the system is decoupled into L-shaped edge pairs. More precisely, S13 7^ 0, sj.4 7^ 0, 
S23 7^ and S23 7^ corresponds to n , j, r and l shaped pairs, respectively, cf. Fig. [T] 
and the boundary conditions ([6]). □ 

Apparently, if Si 2 = S34 = and S contains just one nonzero off-diagonal element 
pair, the spectrum consists of isolated points. 

8.3. Linearly growing spectral gaps, constant bands 

Consider now the case when at least two off-diagonal elements of S are nonzero, 
i.e. the situation when both T>\, £> 2 are intervals. Then the spectral asymptotics is 
determined by (1251) where d is given by ([25]) . Employing the symbols T)i,T> 2 introduced 



High- energy asymptotics of the spectrum of a square lattice 



26 



at the beginning of this section, we may characterize the values k corresponding to the 
spectrum as 

k ^ J — + — + 0(n~ 2 ), — + + 0{n~ 2 ) , 
w \ a nir a mr ] 

3=1,2 \ / 

where dj = minVj, d, = maxVj. Note that d 1 ,d 2 depend only on S, since the term 
(— l) n can be absorbed into 61, 62- The above relation gives spectral values in the form 

k2 e U (^T- + 2 - + ^ + 2 - + > ( 30 ) 

3=1,2 \ / 

which means that spectral bands (or their overlapping pairs) are of asymptotically 

constant width as the band number goes to infinity, and consequently, the spectral gaps 

are linearly growing with the band index. We also remark that the case a nontrivial S' s 

coupling considered in |Ex96, EG96] corresponds to S with all the entries nonzero and 

identical and thus has the described gap behaviour as expected. 

The band structure depends on the parameter values, in particular, the band 
corresponding to T>\ and T>2 in ( 1301) may or may not overlap. The former situation 
occurs, e.g., if only the elements sj.2 and s 34 are nonzero, the latter if Su + s 2 2 — s 33 — s 44 
is large in comparison with the off-diagonal elements. 

Let us finish this section with a short note on the prominent examples of the 5' s 
and 5' coupling. A direct substitution of the appropriate boundary conditions into the 
expressions for dj leads to these results: 

• In the case of the 5' s coupling, it holds d\ = d\ = 0. At the same time, one of 
the sets P 1 ,P 2 equals {d^}. Consequently, the spectrum consists of asymptotically 

A. I 

constant bands and linearly growing gaps between them. We remark that if d 2 = d\ 
and at the same time D\ and T> 2 are non-degenerate intervals, a special analysis is 
needed to find whether there is a gap between the bands corresponding to T>\ and 
T>2 or not. 

• In the case of the 5' coupling, we have d\ < d\, thus the sets T>\ and T>2 are disjoint. 
Consequently, the spectrum has asymptotically the pattern ■ • • GbgbGbgbGbgbG ■ ■ ■ 
where G represents linearly growing gaps and b, g stand for bands and gaps whose 
widths are asymptotically constant. 

Appendix: proof of Proposition 16.21 

We will prove the claims in the order of their presentation: 

(i) It holds -K c cos 2 x + K s sin 2 x = - K -+ K " - Kc + K ° cos 2x and 

max { cos x ■ L c (9i, 82) + L{Q\, #2) | 0i, @2 £ ( — 7T, n]} 

= max{co8x-L c (e 1 ,e 2 ) + L(e 1 ,e 2 )\ 61,62 e (-2tt,0]} 
= max {cos x ■ L c {6\ — ir, 6 2 — tt) + L{6\ — n, 9 2 — it) \ 6\, 9 2 G (— tt, iy]} 

= max{-cosx • L c (6 1 ,6 2 ) + L{6 X ,6 2 ) \ 61,62 G (-tt,ty]} , 



cos2x + max { | cosx| • L c (9i, 9 2 ) + L(9%, 9 2 ) \ 9i, 9 2 £ (— 7r, tt]} , 
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therefore V 2 + (x) equals 

-K c + K S K c + K s 
2 2 

and the function V 2 + (x) can be treated in the same way. To finish the proof of (i), 
it suffices to realize that both cos2x and |cosx| are 7r-per iodic functions satisfying 
cos 2 f| — x) = cos 2x and | cos (| — x) \ = | cosx|. 

(%%) The part of V 2 (x) independent of 9i,9 2 is obviously increasing on [0, |] , 
because — cos2x does. As for the second part, for all x' £ [0, | ) it holds: if 9 1: 9 2 
satisfy 

min {cosx' • L c (9 h 9 2 ) + L(9 h 9 2 ) \ 6 U 9 2 £ (-tt, tt]} = cosx' • L c (9~i, 9 2 ) + L{6 U 9 2 ) , 

then L c (9i, 9 2 ) < 0. If this was not the case one could shift both 9\, 9 2 by tt, which would 
change the sign of L(9\, 9 2 ), and therefore make the value of cosx' ■ L c (9i, 9 2 ) + L(9i, 9 2 ) 
smaller, however, this would contradict the assumed minimality. 

Let now x" < x' £ (0, f ] . Then cosx" • L C {9 U 9 2 ) + L(9\, 9~ 2 ) < cosx' ■ L C {9 X , 9 2 ) + 
L(9\,9 2 ), and therefore 

min{cosx"-L c (# 1 A) + L(# 1 A)| 9 t ,9 2 £ (-7r,7r]} 

< mm {cos x' ■ L C (9 U 9 2 ) + L(9 U 9 2 ) \ 9 X ,9 2 £ (-tt.tt]} , 

thus min-fcosx' • L c (9\,9 2 ) + L(9i,9 2 ) \ 9\,9 2 £ (—71,71]} is an increasing function of x 
in the interval [0, |-] . 

(Hi) Let 9i(x) and 9 2 (x) be functions defined on [0, |] such that 

max {cosx • L c (9i, 9 2 ) + L(9i, 9 2 ) \ 9i,9 2 E (— tt, tt]} 
= cosx ■ L c (0~i(x),0 2 (>)) + L(9 1 (x),9 2 (x)) ; 

it follows that 

8L 8L 
C osx--/(9\{x),9~ 2 {x)) + ^(£(z),0 2 (x)) =0, j = 1,2, (31) 
d9j d9j 

which we will use below. The function V 2 + (x) can be then written as 

— K c + K S K c + K s ~ ~ 

V 2 + (x) = ^ — - cos2x + cosx • L c (9i(x), 9 2 (x)) + L(9i(x), 9 2 (x)) , 

so we can compute its derivative with respect to x, 

= (K c + K s ) sm2x - sinx • LMx), 9 2 (x)) 

dx 

(dL c ~t, . dL c ~i. A dL~i. . 8L ~ 1 . , 
+ cosx — ^9 X x + -^9 2 x + -^0! x + -^9 2 x . 
\d0i d9 2 J d9 x 89 2 



The expression on the second line vanishes due to (l3Tj) . hence 
d 



, -V, fx) = sinx 
dx 2 v 7 



2(K C + cosx - L c (9 x {x),9 2 {x)) 
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Now we will show that < =>■ V 2 + (x) > 1 — once this is proved, the proof 

of statement (Hi) is finished. We have j^V^ix) < =>- — K c cosx > K s cosx — 
^L c (9i(x), 92(x)). A substitution for K c cosx into V 2 + (x) gives the inequality 

V+(x) >K s + -cosx- L c (9\(x), 9 2 {x)) + L(9\(x), 2 (x)) . 

Since K s > 1, it only remains to check that \ cos x ■ L c (9i(x) , 9 2 (x))+L(9i(x), 9 2 {x)) > 0. 
The argument leans on the following two statements: 

51. L c (6 l (x),e 2 (x)) > 0. 

52. If L(§i{x), 9~ 2 {x)) < 0, then cosx • L c (9~x(x), § 2 {x)) > 2\L(9 l (x),9 2 (x))\. 

Statement SI can be demonstrated using the same idea as in part (ii) of this proof, 
S2 can be proved by reductio ad absurdum. Let us suppose that L(8i(x),8 2 (x)) < 
and |cosx • L c (9~i(x), 9~ 2 (x))\ < 2\L(9~ 1 (x), 9 2 {x))\. Then we define 9\ = 9~\{x) + p and 
9 2 = 9 2 {x) + p where p = | or p = — | — the sign of p is chosen such that L c (9i, 9 2 ) > 0. 
Then it holds 

cos x-L c (9\,9 2 )+L(9\J 2 ) > cos x ■ - L(9\(x), 9 2 {x)) = \L(9~x{x),9~ 2 {x))\ 

> ^cosx ■ L c (9\{x),9~ 2 {x)) > ^cosx ■ L c (9 x {x), 9 2 {x)) + L(9\(x), 9~ 2 (x)) , 

hence 

V 2 (x,9 u 9 2 ) > V 2 (xMx)Mx)) = v 2^), 
which is in contradiction with the definition of V 2 + (x). 

(iv) It follows from their construction that the functions V 2 + , V 2 (x) satisfy 
V 2 {x) > V 2 ~(x) for all x G (0, |). The strict inequality V 2 + (x) > V 2 (x) is equivalent to 
the fact that the image of the function cosx • L c (9i, 9 2 ) + L(9\, 9 2 ) defined on (— n, 7r] 2 
forms a non-degenerate interval. Since cosx • L c (9i, 9 2 ) + L(9\, 9 2 ) equals 

9? (-8 cos ak(t n t^ + t 12 t^)e 1&1 ) + 3? (8 cos ak(t 22 t^ + t^t 12 )e W2 ) 

i.e. it is of the type examined in Proposition 16.11 we infer that the image degenerates 
to a single point iff 

^11^21+^12^22 = A £ 2 2*2i + *ii*i2 = A t u t^ = A t 12 t^ = ; 

it is straightforward to check that these four conditions are satisfied iff at most one of 
the numbers tn,ti 2 ,t 2 i,t 22 is nonzero. This concludes the proof. 
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